RECONSTRUCTION OF LESS REGULAR CONDUCTIVITIES IN 

THE PLANE 
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Abstract. We study the inverse conductivity problem of how to reconstruct an 
isotropic electrical conductivity distribution 7 in an object from static electrical 
measurements on the boundary of the object. We give an exact reconstruction 
algorithm for the conductivity 7 G C^~^'^{^1) in the plane domain il from the asso- 
ciated Dirichlet to Neumann map on dVl. Hence we improve earlier reconstruction 
results. The method used relies on a well-known reduction to a first order system, 
for which the 9-method of inverse scattering theory can be applied. 



1. Introduction 

Let Q C be a bounded and smooth domain and let 7 G L°°{fl) with < c < 
7(0;), a.e. X E Q for some constant c. Consider the conductivity problem 

V ■ 7VM = 0, in n, 

(1) 

M = /, on dQ, 

for / G H^^'^{dVL). This problem has a unique solution u G H^{VL) and hence we can 
define the Dirichlet to Neumann map 



an, 



where u is the outer unit normal defined in a neighborhood of dQ. This definition 
should be understood in the sense, that if f,g E H^^'^{dfl), then 



(A7/, g) = 7VM ■ Vvdx, 
Jn 



where u G H^{Q) solves (P and v G H^{Q) has v\gQ = g. 

In this paper we study the problem of reconstructing 7 from knowledge of A^. This 
is the reconstruction issue of the inverse conductivity problem. 

The main motivation for studying this problem comes from a method for medical 
imaging called Electrical Impedance Tomography (EIT) (see ||CIN99|| for a thorough 



review of current methods in EIT). In EIT one has a conductive body with un- 
known conductivity 7, and from static electric measurements on the boundary, i.e. 
by applying a voltage potential on dQ and measuring the current flux through dfl, 
one would like to monitor the interior conductivity. Since tissue, bones, lungs etc. 
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have different conductive properties, an image of the conductivity is an image of the 
internal structures. 

The mathematical formulation of the inverse conductivity problem was introduced 
by Calderon in [pal80|| . Calderon gave an approximate reconstruction method for 7 



close to constant. Since Calderon's fundamental paper a lot of work has been done 
on this problem. 

The first global result was in ||KV84|| , where it was proved, that if 7 G C°^(fi) then 
determines 7|an and all partial derivatives on dQ. From here it follows by analytic 
continuation, that if 7 is real-analytic in Q then it is determined by A-,. Later in 



KV85II the result was extended to include piecewise real-analytic conductivities. 



In three dimensions a major breakthrough in the theory of inverse boundary 
value problems was the paper ||SU87|| . There it was proved that A^ determines 
7 G C°°{Q). Their method of proof relies on the existence of special solutions to 
the conductivity equation, the so-called complex geometrical optics solutions or ex- 
ponentially growing solutions, which behave like e^'^ for large |x| and for C G C" 
with ( ■ ( = 0. The assumption on 7 has since been relaxed by a number of au- 
thors ( [|Ale ^ 



| Bro96|| ). The best uniqueness result so far is for 



7 



G obtained in | [FPU[| . Using exponentially growing solutions Nachman 

( [[Nac88|| ) gave a reconstruction algorithm for conductivities 7 G C^'^(n). Further- 
more, Alessandrini ( |[Ale88|| ) proved a conditional stability result. 

In two dimensions, uniqueness for small 7 G W^'°°{Q) was proved in ||SU86|| , but the 
general question remained open until 1996, when Nachman |[Nac96|| gave a uniqueness 
proof and a reconstruction algorithm for 7 G W'^'''^{Q),p > 1. Nachman's method 
relies on the reduction of the conductivity equation to a Schrodinger equation, which 
essentially requires two derivatives of the conductivity. For this equation he proved 
existence of exponentially growing solution, which can be reconstructed from A^ by 
the so-called 5-method in inverse scattering theory originally developed by Beals 
and Coifman (see [ [BC89|| and the reference there). Then using the exponentially 
growing solutions Nachman was able to reconstruct 7. The algorithm has been tried 
out numerically (see |[SMlUCm ). In 1997 Brown and Uhlmann (| |BU97| |) improved the 
uniqueness result for 7 G W^'^{Q),p > 2. Their method relies on a reduction to a first 
order system for which only one derivative is required. For this system the (9-method 
is again applicable, see Section |^ for further details. For both works conditional 
stability has been proved, see ||Liu97|| and ||BBR01| . 

For a thorough review on the history and developments of the inverse conductivity 
problem we refer to ||Uhl99 . 



In this paper we will show how the uniqueness proof in |[BU97|| for the two- 
dimensional inverse conductivity problem can be turned into a reconstruction method. 
The reconstruction algorithm parallels Nachman's algorithm, but since the reduction 
used here requires less regularity, we can improve the reconstruction result in ||N ac96 . 
The main result is the following: 

Theorem 1.1. Let Q G M."^ be a bounded and smooth domain. Let e > and let 7 G 
C^~^''{fl) be such that < c < ■y{x) for some constant c. Then 7 can be reconstructed 
from Aj. 
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The proof is given in the following sections as follows: In Section |^ we introduce 
notation, outline the main idea of the uniqueness proof and give the reconstruction 
algorithm. The steps of the algorithm are then justified in the following sections. 

2. Outline of the reconstruction 

In this section we introduce notations and review the reduction of (JlD to a first 
order system. We then formulate an inverse problem for this system and finally 
outline the reconstruction method. 



Generally we work under the hypothesis stated in Theorem |1 . 1| . Furthermore, we 



assume that 7=1 near dil. It was noticed by Nachman in |Nac96|, that this assump- 



tion is not a restriction. When 7 G C^~^'^{Q) it is well-known ( [ |Ale8 8|| , | |S U 8 8| | , | [N ac8 8|| ) 
that 7 and its normal derivative on the boundary are determined by A-y. Hence 7 
can be extended outside Q, such that it is identically 1 outside a compact set and 
7 G C^"'"'^(R^). Moreover, the Dirichlet to Neumann map on the boundary of a larger 
domain containing the original domain can be reconstructed from the original Dirich- 
let to Neumann map and from the known extension. We show in the appendix how 
to do such a construction for 7 G L°°{Q). 

It is well-known from the theory of elliptic equations, that if / G C'^~^'^{dQ) and 
7 G C^^^i^fl), then (|1|) has a unique solution u G C'^'^''{fl). Then the vector valued 
function {v,w)'^ defined by 



du 



(2) g=y/=^^;j.c-(n)xc-(n) 

satisfies the equation 

(3) (D - Q) = 0. 

Here the matrix operator D and the matrix potential Q are defined by 

where d = {d^ — idy)/2 and d = {d^ + idy)/2 and 
(5) q = -7-V25^i/2_ 



We shall see below in the proof of Theorem |3]^ that there is a one-to-one corre- 
spondence between u and (f , if )^, i.e. given a solution {v,w) to (H), where Q is of 
the form with q given by (|^), there exists a solution u to (|l]) such that (0) holds. 

With Q G C"^(]R^), the natural Cauchy data for the system is 

Cq = {(i^loc^loc) : {v,w) G C^+^p x C^+^p, {D - Q){v,wf = 0}. 

We can now pose the following inverse problem for the system (^): Let q G C'^{VL) 
and Q defined by (|). Is Q uniquely determined from Cq and, in case it is, how can 
it be reconstructed? 

In this paper we show how to reconstruct Q coming from a conductivity equation. 
For this particular case we give an explicit characterization of the Cauchy data Cq in 
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terms of the Dirichlet to Neumann map and this characterization enables us to 
solve the inverse problem. 

Brown and Uhlmann ( ||BU97]| ) were the first to use in the study of the inverse 
conductivity problem and their uniqueness proof was based on the scattering theory 
for this system developed by Beals and Coifman ( ||BC88|| ). For the direct scattering, 
one looks for a family of functions, which solve the equation (J^) in the whole plane 
(with Q extended by zero outside f2) and have a certain exponential behaviour at 
infinity. More precisely, for each G C, one seek solutions '^{■,k) of the form 

(6) ^(z, k) = m{z, k) 









~izk 



where m approaches the identity matrix / as \z\ ^ oo and z = Xi + 1x2 G C. These 
solutions are called exponentially growing solutions or complex geometrical optics. 
A simple calculations shows that m should solve the equation 

(7) {Dk - Q)m = 0, 

where Dk is the matrix operator 

'(d - ik) 
{d + ik)^ 

The analysis of the direct scattering problem along these lines is due to Beals and 
Coifman ( ||BC85| |, pUg^ ). A detailed analysis for the case Q G ^^(M^) n L°°(M2) can 
be found in |{Sun94a| | , | [Sun94b| | , | |Sun94c| ] , and the analysis for the case of compactly 
supported Q G LP(M^),p > 2, can be found in jBU9ti|]. G iven our assumptions on Q 
we have the following theorem, see pU9q| and IPBROlj . 

Theorem 2.1. Let Q he a hermitian, compactly supported matrix with Q G C"^(M^). 
Then for any 2 < r < 00 and k E 'C there is a unique solution m{z, k) to (|7|) with 
(m(-, k)-I)e L"(M2) n C'iR^). 

Furthermore, for fixed z G M^, the map k \—>- m{z, k) is differentiable, and it satisfies 
the equation 



— m(z, k) 
dk ^ ' 



m{z,k)Ak{z)S{k), 



where 

(9) 

and 



S{k) 





-q{z)muiz, k) 



q{z)m22{z, k) 




dfi{z) 



e{z, k) 






e{z, —k) 



for e{z, k) = exp{izk + izk). 



In the formula mn and m22 are the diagonal components of m, and dfi{z) is 
the usual area measure. Note, that since Q G ^(M^), elliptic regularity theory gives, 
that m G C/+^(M2)_ This fact will be used later. 
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The matrix S{k) is called the scattering matrix. Using the definition of \E' and the 
fact that Q is supported in Q, an integration by parts in (Rl) shows that 



(10) sw = f / ( .3- °, , » "(f 

where il)i2 and ip2i are the off diagonal entries of In the above formula, da{z) is the 
Lebesgue measure on the boundary and u is the complex normal at the boundary, i.e., 
if (z^i(^), i^2{z)) denotes the outer unit normal at z G dfl, then u^z) = I'li^z) + ii'2{.z) 
and V = I'llz) — iv2{z)- Therefore, in order to determine the scattering matrix S'(fc), 
it suffices to find the traces on dVl of the exponentially growing solutions. 

Once S{k) is known, we can solve the 9-equation (H) to find m(z, k) for each fixed 
^ G C and reconstruct Q by the formula 

Q{z)= lim /i(5,(0))"i / Dkm{z,k)dfi{k) 

^0-*°° J {k:\k-ko\<r} 

from ||BU97|| . This would solve the inverse problem for (^). Furthermore, 7 can be 



found by solving (^. 

Since we are mainly interested in reconstructing the conductivity, and this approach 
then requires taking the parameter to infinity and solving a boundary value problem, 
we shall take a different route. As observed in | |BBK01|| we can reconstruct 7 directly 
from m{z, 0). In fact we will see, that 

(11) I'^^z) = Re(mii(z, 0) + m2i{z, 0)), 



where now m is the unique solution given by Theorem to Dkm = —Q'^m. It was 
noticed by Reals and Coifman that the scattering matrix S associated the potential 
—Q^ relates to the scattering matrix 5* for Q by the formula 

(12) S{k) = S{-kf- 

Therefore rh can be found from S by solving (P) with 5* replaced by S{—k)'^. 
We are now ready to outline the reconstruction method: 

1. Determine \E' on dQ (see Theorem |4.3| below) 

2. Define the scattering matrix S by (p!OD and S by (|T2|) 

3. Solve the 9-equation (H) with S for m{z, /c), z G 

4. Recover 7 on i7 from (|n|). 

The implementation of this method is subject to further study. 

3. A BOUNDARY RELATION 

In this section we give an explicit characterization of the Cauchy data for the first 
order system (H), in case the potential Q comes from a conductivity. 

Assume u G C^+^U) solves (g) for some / G C^+'{dn). Then {v,w) = -f^/^{du,du) 
solves (Bf) and one easily finds the relation 



(13) 



on dQ, 
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where ds denotes the tangential derivative on dQ. Inverting (|T3|) gives 



(14) 



an 



9s{f) J \iv -iv) 

Let : C^^'^{dVL) —>■ C'^^'^{dQ) denote the normalized inverse of dg on the boundary 
defined by 

d;'f{s)= r fds + c{f), 



where the measure ds is the arc length on dfl and the constant c(/) is chosen so that 
Jg^d~^ fds = 0. Then from (|l^) we arrive at the following boundary relation for 

This relation motivates the following definition 

Definition 3.1. Let (hi, he a pair of functions in C^^''{dVL) x C^^'^{dVL). We say 
that (/ii,/i2) satisfies the boundary relation if 

(15) iK^d'^^uhi — Vh2) = {vhi +Vh2). 

We just saw that solutions to (||), defined through (||), satisfy the boundary relation, 
but as the following theorem shows, (0) is in fact a complete characterization of Cq : 

Theorem 3.2. If Q is given by (|^) then 

(16) Cq = {{hi, h2) e C^+'{dQ) X C^+'{dQ) : iA^d:\uhi - vhi) = uhi + Vh2}. 

Proof. First we show that any pair {hi, h2) G C^^^{dVL) x C^^^{dVL) which satisfies the 
boundary relation is in Cq. Let u G C^^^(f2) be the unique solution of the Dirichlet 
problem 

'V-7Vu = 0, inn, 

u = id'g^ {uhi — Vh2) , m. dVt. 

Define a solution {v, w)^ to (H) by the relation (0) with u from above. That {v, w) = 
{hi, h2) follows from (^) by the use of ([151) . 

Conversely, we start with a solution {v,w) G C^^'^{VL) x C^^^{Vt) to (|^) and check 
the boundary relation for {v,w)\qq, . Since v and w solve the linear system (^ and 
Q is of the form (^), we have the compatibility relation 

%-V2^) = d{-i-^'^w). 

The Poincare lemma ensures now the existence of such a u that 



du 



w I \du 



It is easy to check that -u is a solution to the conductivity equation in the form 
2dud'~^ + 2(9m(97 + i'jddu = 0. Now relation (|13D with / = u\q^ shows that {hi, h2) 
solves the boundary relation. □ 
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4. From to S{k) 

In this section we show how to reconstruct the trace on dQ of the exponentially 
growing solutions "^{z, k) defined by (|^). Notice that the first row has entries which 
are analytic outside Q, while the second row has entries which are anti-analytic 
outside Q. Moreover they have a prescribed behavior at infinity and their traces 
on dQ satisfy the boundary relation (|TB|). We will prove, that these relations are 
sufficient to reconstruct the trace of uniquely. 

Due to the symmetries 



(17) mii{z,k) = m22{z,k), m2i{z, k) = m^iz, k), 

which follows from the differential equations, the asymptotic for the columns of m and 
the uniqueness in Theorem it suffices to reconstruct the first column (^'n, ^1^21)^ 
of k). 

For every {(, z) E C with ( ^ z we introduce 

I g-«fc(C-2) 

9kiC,z) = , 

IT Q — z 

a Green's kernel for d which also takes into account exponential grow at infinity. 
Using (yffc(-, z) with z G dVt, we define the single layer potentials Sk and Sk as boundary 
integral operators by 

(18) S,f{z)= [ f iOg.iC, z)dC Skf{z)= [ f{C)UC.z)dC. 

These integrals are understood in the sense of principal value. It is a classical result 
in singular integral theory that these operators are well defined, see ||Mus53|| . 

Proposition 4.1. //O < e < 1 then Sk, Sk : C^^'{dn) C^+'{dn) are well defined 
bounded operators. 

Let /C denote the matrix operator 

1 fl-iSk 



^ 2V I-Sk. 

The following lemma gives a necessary condition for a pair of functions defined in 
Qf^ = C\fi to be analytic and anti-analytic respectively and have a special exponential 
decay at infinity: 

Proposition 4.2. Let v he analytic and w be anti-analytic in Q'^ and continuous to 
dVl, such that ve^''^^ — 1 G U{VL'^) and we~^^'' G for some r > 1. Then the 

trace (/ii,/i2) = {v,w)\gQ solves 

(19) /C'^^ ' 



h2j \0 

Proof. We will prove the relation for v; for w the result follows by similar arguments. 
Let z E Q'^ he a fixed point inside some large ball B^i of radius R. With the d 
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operator acting with respect to C and R > such that Q C B{0,R), Pompeiu's 
formula ( [|Vek62|] ) states, that 



v(z)e-"'' - 1 = / dC / dC 

d{v{C)e-'^^ - 1) 

at, 



27ri Jqbr C- z "^T^i Jdn C - ^ 

Since fe~*^'^ — 1 G L''(]R^) the integral on dBfi converges to zero at least on a sequence 
of increasing radii, and since z G Q'^, Cauchy's formula gives 

1 f 1 



dC = 0. 



27rz JanC- z 
Therefore 

(20) v{z) + i / viOgkiC, z)dC = e^^^ z e Q''. 



2i 



an 



Now let z & Q'^ approach some point on the boundary and apply Plemelj's Formula 



Mus53| to get the equation for f . □ 



We know that the exponentially growing solutions {ipii{-, k),ip2i{-, k))\an must sat- 
isfy (19) and the boundary relation ([T5|). The following result shows that these 
conditions are sufficient to characterize the trace of the exponentially growing solu- 
tions. 

Theorem 4.3. The only pair {hi,h2)^ G C^^''{d^l) x C^^''{d^) which satisfies 
(21) 



ho \ 



(/ - ik^d:^){vhi){z) = (/ + iK^d;^){Vh2){z), 
is (V^ii(-,/c),V^2i(-,fc))|an- 

Proof. Let h = (/ii,/i2) G C^+'{dn) x C^+'{dn) be any solution to (gT]). Guided by 
and the related formula for w we extend h outside Q to {v{z), w{z)) by 

v{z) = -^[ h,{09k{C,z)dC + e''\ yj(,) = -Lf h2{0UC,z)dC. 
2« Jan 2« Jan 



Inside fl we extend /i as a solution to (0), which can be done by Theorem since 
h satisfies the boundary relation. 

We will prove, that {v, w) is a solution to (|^) and that ve'"^^^ — !, w G L''(M^), r > 2. 
Then the result follows from Theorem |2.1| . 

Outside Q, v is analytic and w is anti-analytic. Moreover, since 7 is identically 
1 near the boundary, we also have that v is analytic and w is anti-analytic in an 
inside neighborhood of the boundary. We need to check the analyticity of v and anti- 
analyticity of w as we go across the boundary. For this it suffices to have continuity 
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across the boundary and then invoke Morera's theorem. We shall show the continuity 
of f , for w similar reasoning works. Let z approach some point zq G dVL from outside. 
Then using Plemelj's formula we get 

\^^v{z) = - (-^ + ^^S,h{zo)) + e-o^ 



Now use the first equation of (|2TD to conclude \imz^zQv{z) = hi{zo)- The continuity 
of V from inside comes from the regularity theory for elliptic equations. 
We are left to show that fe~*^*^ — 1 G L^(M^), but since 



v{z)e-^^' - 1 = -— / ^ dC = O — , as 1^1 oo, 

it is in U{E?) for r > 2. This proves the theorem. □ 

As an corollary we easily get 
Corollary 4.4. ker /C n Cq = {0}. 

Proof. Let {hi, h2)^ G ker /CflCg and (f , wY be the trace of the exponentially growing 
solutions. Then {hi + v,h2 + w)^ is also satisfies the boundary relation and satisfies 
(|2T|) . Now use the uniqueness of Theorem ^]3| to conclude that {hi, h2)^ = 0. □ 

To obtain the exponentially growing solutions on dQ one has to solve the three 
equations (^TJ). Since {e'^^^,0)'^ solves 



/C 



\ 



the pair of traces of the exponentially growing solutions can be constructed as the 
projection of {e^''^,0)^ on Cq along ker/C. 

5. From S{k) to 7 

In this section we will prove how to reconstruct 7 from S{k). More precisely we 
will work with the dual S{k) = S{—k)'^, solve the equation @ for m and finally find 
7 by the formula (pJ]). 

Since the equation (|]) is mixing terms, we will start by writing pseudoanalytic 
equation for a suitable linear combination of the elements in m : 

Lemma 5.1. Let z EC be fixed and let 



(22) m+{z,k) = mii{z,k) + mi2{z,k). 

Then for k E C we have the equation 

d ^ 



(23) -^m+{z, k) = e{z, -k)S2i{k)m+{z, k), 

Proof. The result follows from (^ and (0). □ 
To solve this equation we will need a result about the behavior of m+ for large k : 
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Theorem 5.2 ( ||BU97|| ). Suppose Q G LP{M?),p > 2 is compactly supported. Then 
for all q > 2p/ {p — 2) 

sup \\m{z, ■) - /||l<7(k2) < C 

zee 

where C depends on q and Q. 

Furthermore, we wiU need a certain decay of S. hi ||BU97|| the fact that S G L^(R^) 
was sufficient to prove uniqueness of solutions to ( p3D subject to the asymptotic given 
in Theorem |5.2| . We will however, need a slightly better decay for 5* in order to write 
the integral equations corresponding to (pS]). The following lemma from | BBR01 
gives the result: 

Lemma 5.3. Let Q G C^(R^) he compactly supported. Then S G L''(R^) for any 
r>ro = 4/(2 + e). 

To simplify notations we introduce the operator {d/dk)~^ = dj^ defined by 

We can now write the integral equations for m+ and prove unique solvability of these: 

Lemma 5.4. Let Q G C"'(M^) be compactly supported, and let z & C be fixed. Then 
for g > 4/e the equations 

(24) (/ - d^\e{z, -k)S2r)){m^ - I) = d^\e{z, -k)S2i), 

have the unique solution m+{z, ■) — I in L'^(M^). Furthermore, mj^{z, ■) — / G C"(M^), 
fora< (1 +e)/2. 



Proof. Since £ L^(]R^) and \e{z,—k)\ = 1 we know (see ||Nac93| , Lemma 4.2]) 
that d^^{e{z,~k)S2v) is a compact operator in L''(]R^) for 2 < s < oo. Furthermore, 



since S21 G L''(R^) for some 4/(2 + e)<r<2 (Lemma |5.3| ), it follows by the Hardy- 
Littlewood-Sobolev inequality ( ||Ste7q| ) that d^\e{z, ~k)S 21) G ^^(R^) for q > e/2. 
Hence the integral equation ( P^ in L'^(R^) follows from ( PBj ) and Lemma ^72 . 

The unique solvability of {U4) in -^v'^(R^) follows from Fredholm alternative if 
(/ + d^^{e{z, —k)S2i-)) has trivial kernel in L'?(R^), but this is known from Liou- 

ville's theorem for pseudoanalytic functions since e{z, —k)S2i G L''(R^) nL'"'(R^) (see 
I BU91 , [[Sun94a|| , ||N ac96|| , [|Vek62|| ) . 

To prove the Holder continuity of rhj^{z, •) — / we use the fact, that convolution by 
1/z maps LP(R^) into C"(R^) for 1 < p < 00, a = 1 — 1/pq and max(2,p) < po < 00, 



see 



By the previous theorem we find for fixed 2; G C 

(m+(;., fc) -/) = (/ - d^\e{z, -k)S2v))-\d^\e{z, -k)S2i)), 
in L'^(R^) n C"(R^) for a = 1 — (1/r + l/g). To reconstruct 7 we use the formula 



7(z) = mii(z, 0) + m2i{z, 0) 
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from ||BBR01| , Proposition 4.2]. Since 7 is real, we rewrite the formula using (|T^ and 



(22) 



7(z) = mu{z, 0) + m2i{z, 0) 



This ends the reconstruction. 



Re mil {z, 0) + Re m2i {z, 0) 
Re mil (-2)0) + Re mi2 {z, 0) 
Rem_|_(z, 0). 

Appendix 



We will prove, that the Dirichlet to Neumann map for an extended domain Q' can 
be calculated from the Dirichlet to Neumann map for the domain Q given the value 
of the conductivity in fi' \ : 

Lemma A.l. Let 7 E L°°(fi) and let be given. Let Vl' C be such that C ^l'. 
Take any function 7' G L°°{Q') with 7'ln = 7- 
Then Ay can be recovered from A^ and 7'|f7'\n 

Proof. From the definition of the Dirichlet to Neumann map we have for any f,gE 
H^/\dn'), that 

(Ay/,(?)= [ yVu-Vv 
Jn' 

= _'y'Vu -Vv + {A^u\dn,v\an), 
Jn'\n 

where u G H^{fl') denotes the unique solution to 

V ■ 7'Vm = 0, for X G Q', 
u = f, for x G dQ' , 

and V G H^{Q') is any function with v\dQ' = g. Hence we see, that Ay can be 
calculated from 7!^/^^ ^7 without exphcit knowledge of 7 in fl, provided the 
solution u to (^) can be found in Q' \ Q. 

We claim, that m in fi' \ i7 can be found as the unique solution to 

V-yVM = 0, fora;Gfi'\n, 
^2Q-j u = f, for X G dQ', 

Ou 

^'o^ldo = A^{u\on), for x G dn' . 
That mIq/^q solves (^) follows, since for any v G Hq{Q') 

/ _'y'Vu -Vv + {A^u\an,v\an) = / 7'Vm ■ Vf 
Jn'\n Jn' 



0. 
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For the uniqueness we assume, that uq G H^{Q' \ fl) solves with / = 0. Now 
extend uq into Q as the solution to 

V ■ 7VM = 0, for X eQ, 

u = uo, for X G dQ. 



But since Uq is in H^{Q') and solves 



with / = 0, we conclude, that Uq 



in 
□ 
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